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Abstract 

Real time excitations in the broken symmetry phase of the classical Abelian Higgs model are 
investigated numerically in the unitary gauge. Spectral equations of state of its constituent 
quasi-particles are extracted. Characteristic differences between the statistical correlations of 
the transverse and of the longitudinal vector modes with the Higgs field are exposed. A method 
is proposed for the reconstruction of the quasi-free composite field coordinates. 

1 Introduction 

The rearrangement of the multiplet structure of elementary fields in Gauge+Higgs models in the pro- 
cess of symmetry breaking is one of the most fundamental concepts upon which the Standard Model 
of elementary particle interactions is built PlEj- The transmutation of the pseudo-Goldstone field 
into the longitudinal polarisation state of massive vector bosons at the minimum of the effective po- 
tential has been successfully used in perturbative calculations. The absence of light (massless) fields 
from the spectra of the electroweak sector of the Standard Model was also tested non-perturbatively 
first at relatively strong couphng [3], later close to the continuum lUE]. Real time investigations of 
classical gauge field dynamics mostly concentrated on topological aspects of the process of baryogen- 
esis, e.g. measuring Chern-Simons densities jElIZI, and on the emergence and the evolution of gauged 
strings [Hlinilini- 

The idea of far-from-equilibrium baryogenesis [TT| IT^ IT^ ITU IT^j led to increasingly detailed 
studies of the real time mechanism of the reheating of the Universe at the end of inflation. Features 
of the parametric resonance and tachyonic instability were first studied in simplified infiaton-|-Higgs 
systems fU ^ ITT f IH^ j realising the hybrid inflation scenario. In these systems the breaking of the 
global symmetry leads to the generation of Goldstone fields and related global topological excitations 
(e.g. strings). Recent numerical studies of Higgs systems started the detailed quantitative exploration 
of the real time excitation process of the Higgs and massive vector particles |20[ I15| ITUj. 

J. Smit et al. jSOl EHj developed methods for the identification of the different physical degrees 
of freedom in the process of excitation in Coulomb and unitary gauges with help of measuring the 
corresponding dispersion relations. They measured the time variation of the corresponding number 
densities during tachyonic (spinodal) instability. In the unitary gauge they find significant differences 
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Dexween tne eariy values xaKen oy ine numoer aensiiies oi lue transversal ana longixuamai aegrees 
of freedom. The differences diminish when the system approaches equihbrium. 

Our final goal in the real time investigation of the Abelian Gauge+Higgs system is to measure 
the evolution of the equations of state of its physical components during the reheating period. To 
achieve this we apply methods developed in our earher investigation of pure scalar (inflaton+Higgs) 
systems |211 119j . In order to understand the nature of the final state to which the system eventually 
relaxes in this paper a careful numerical study of the real time fluctuations in equilibrium systems 
with low energy densities is presented. We extract the equilibrium equations of state of the differ- 
ent quasi-particle species and an attempt will be described to establish the corresponding effective 
field-coordinates. Such dynamical variables are present in rather strongly interacting systems of con- 
densed matter physics, even when the original particles might decay into each other. Of course, this 
investigation could be performed also in the Euclidean formulation of equilibrium field theories on 
lattice. For the clarification of specific issues already at this stage the possibility of real time studies 
of the non-equilibrium time evolution was quite useful (see below). 

We solved numerically the field equations of the Abelian Higgs model, discretized on a spatial 
lattice in the axial {Aq = 0) gauge. The initial conditions were chosen with very low energy density 
in order to end up in the broken symmetry phase. In the present investigation mainly the couplings 

= 1, A = 6 were used. The influence of the variation of the couplings on our findings was 
systematically explored in a broad range. 

As a simple (global) signature of reaching near thermal equilibrium we have chosen the equality of 
the energy density stored in the scalar field (defined as part of the Hamiltonian density independent 
of the gauge fields) with one third of the energy density in gauge modes. It took surprisingly long 
time to reach stationary equality: only after t 4 x 10^|m|~^ did fluctuations in the difference of 
the energy densities settle at the one percent level. Then characteristic thermodynamical quantities 
were analysed for an extended time interval, sufficiently long to explore the thermal features. It is 
remarkable that the equations of state, to be discussed below, take their equilibrium form considerably 
before the true equilibrium is reached. This experience supports the findings of |^ ■ 

The equilibrium configurations were transformed to the unitary gauge and decomposed into a 
scalar (Higgs) field, plus transversal and longitudinal gauge fields. The full energy density and 
pressure of the system can be split up in a rather natural way into pieces associated with these fields. 
In section 2 we shall present the method of constructing spectral decompositions for these quantities 
without making any assumption about the nature of the statistically independent field variables. 
Still, we can convincingly argue that the quasi-particle interpretation of the thermodynamics works 
and the extracted mass degeneracy agrees with the result of the weak coupling analysis. 

In section 3, evidence will be shown that the transverse vector modes behave as independent 
Gaussian statistical variables. On the other hand, the longitudinal vector mode, defined with the 
usual projection violates the expected Gaussian relation between the quadratic and quartic moments 
and does not average independently from the Higgs mode even when the average energy per degree 
of freedom was one thousendth of the Higgs mass. The possible dependence of this characteristic 
difference on the Higgs vs. vector mass relation (that is on the couplings e^. A) was carefully explored. 
Also the effect arising from the presence of vortex-antivortex pairs was understood with help of 
studying the non-equilibrium time evolution of the system. 

The persistent interaction between the longitudinal gauge and the Higgs field motivates the search 
for the quasi-particle "coordinates" in terms of non- linear combinations of these fields. Section 4 
presents the results of the first step made in this direction. In our search for a Gaussian field related 
to the longitudinal vector mode, and fluctuating independently from the other three fields, a new 
pair of canonical variables is proposed: 

Ai\^,t) = (1 + V^x,t))Ai(x,t), PP(x,t) = \ ,.. n^(x,t). (1) 

1 + ap^ x,t 
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ine vame or trie coemcrent a was actcrmrnea oy tne requrrement to nave ror J\l\I l) ana p rnaepen- 
dent statistics. Having done this the new variables are checked for Gaussianity. Also an alternative 
composite field of the form 

Af{^,t) = p(x,t)"Ai(x,t), (x,t) = p(x,t)-"ni(x,t) (2) 

will be tested with satisfactory results. The change of variables influences only very mildly the 
Higgs fleld itself. This is rather fortunate since this variable, as well as the transverse gauge degrees 
of freedom, consistently obeys free quasi-particle thermodynamics. Further possible tests of the 
existence of a Gaussian composite field involving the longitudinal vector mode are discussed in 
the Conclusions. Technical details of the real time numerical investigation are summarized in an 
Appendix. 



2 The thermodynamics of the Abehan Higgs model 

In the first part of this section we review the continuum expressions of the relevant thermodynamical 
quantities of the Abelian Higgs model in the unitary gauge. A model-independent spectral test is 
performed which confirms that in the time interval where the analysis of the thermal features is 
realized the energy is equally partitioned in momentum space among the different modes. After 
checking that the potential energy density and the kinetic energy of the spectral modes are equal in 
equilibrium this fact is used to derive simple expressions for the spectral equations of state determin- 
ing the thermodynamics of the Higgs, the longitudinal and the transverse vector modes. The mass 
degeneracy of the different vector polarisations is demonstrated. The corresponding effective masses 
are close to the tree level expectations. 

The Lagrangian of the Abelian Higgs model is given by the expression 

L = -\f,,F^''' + - V{^), (3) 

where F^j, — d\p^Aj/^, 0^ — 8^ + ieA/j,. V{^) is the usual quartic potential of the complex Higgs field 

= Im2|$|2 + A|$|4^ m2<0. (4) 

The result of the standard calculation of the energy-momentum tensor for the Abelian Higgs model 
leads in the unitary gauge to the following decomposition for the energy density: 



e — Sp + St + cl, 
2^2 



1, 



ST = ^[n^ + (V X At)' + eV'A'], 

e.-^[ni + eV(Ai + A^)]. (5) 

Here the indices T (L) refer to the fact that the corresponding terms contain the transverse (longitu- 
dinal) component of the vector field and its canonically conjugated momentum. The energy density 
is written in Hamiltonian formulation employing the relations: 

Up^p, Ut^At, Ul^Al + VAo. (6) 

Also one has to understand as a dependent variable in view of the Gauss-condition: 

Vn^ = e^p'Ao. (7) 
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ine space-space aiagonai componeras oi ine energy- momemum tensor give simuar expressions lor 
the partial pressures: 



P = Pp + Pt+Pl, 

Pp = ^nJ-i(Vpf -np), 

p^=l[n| + (Vx A^)2-eVA|], 
o 

p^ = ^[ni-eVAi] + l-l,(vn,)l (8) 

Spectral energy densities can be introduced for each of the three contributions defined in Eq.(jHI) 
by considering the square root of each ei(x, t), i = p,T,L and taking the absolute square of their 
Fourier-transform: 



e.M= i4\,e^Wl"\Kt). (9) 



(2vr) 



There is some ambiguity in factoring the densities into two identical powers. This choice is the 
most natural for free quasi-particle interpretation of the thermodynamics, since those expressions 
are quadratic in the effective field-coordinates. In the broken symmetry phase at equilibrium one 
expects |ej-^^^''(k, t)p to fluctuate around k-independent values obeying the relation: 



|ei'/'^(k)|2 : |ei.'/'^(k)|2 : \e^p'^\\i)\^ = 1:2:1. (10) 

The overline means an average over the different modes with the same length fc, and over a certain 
time-interval. This was found to be satisfied with 10% relative fluctuation in the relevant time 
interval as illustrated in Fig. ^ When below we refer to a sort of "temperature" we have in mind the 
average energy density in the k-space. If a free quasi-particle model works, then this temperature is 
twice the average kinetic energy per mode. 

Similar construction can be performed also for the partial pressure of the z-th constituent, pro- 
viding a model-independent definition for its spectral decomposition: 



p-(k)^bf/^)(k)p. (11) 
The central step in the analysis is the investigation of the spectral equation of state defined as 

..(k) . ||. (12) 

This quantity can be measured for each configuration. In the equilibrium the spectral densities of 
the kinetic and of the potential energy are equal on average. Using this feature as a hypothesis for 
the Fourier-modes a very simple functional form is predicted in case, for instance, of the Higgs-mode: 



^ 1 k'|p(k)p 

3k2|p(k)P + 2V,(k)' 

If the Higgs-field is a "pure quasi-particle coordinate", then it performs small amplitude oscillations 
near its equilibrium value p^: 

V{p) = \Ml,^f{p-p,)\ V,{\.) = \Ml,ff\p{\^)\\ MO. (14) 
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Figure 1: Energy densities of the Higgs and gauge fields in equifibrium, corresponding to etotai — 
3.77|m|^. The data points and the error bars were obtained by binning and forming the mean energy 
per mode and its fiuctuation for each bin. 

This leads to a one-parameter "free- field" form of its spectral equation of state: 

(15) 
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Similar line of thought can be followed also in the spectral analysis of the transversal gauge energy 
density and pressure. The result is the following: 



If one finds a good description of the measured spectral equation of state by this formula with a k- 
indcpcndent mass, this provides evidence for the free quasi-particle interpretation of the transverse 
part of the thermodynamics. 

For the longitudinal polarisation a very similar formula can be obtained when the ratio ^^(k) is 
defined by the scaled energy density e^p^eL and pressure e^p^pL as 



tWr k = ^ 17 

' |(epv/iZ)(k)|2 ^ > 

It leads to the above quasi-particle form if the mode-by-mode equality of the scaled kinetic and 
potential energies is obeyed in the form 



k^|nak)|2 + |(epn^)(k)P = |(eVAL)(k)|2, (18) 
with the following effective mass formula: 



Hrfk) 
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Figure 2: Spectral equations of state computed from Eqs. (fT^ and (fTTj) as applied to the three 
constituting terms of the energy density and pressure. 

The first task was therefore to check numerically if the average equality of the kinetic and potential 
mode energies is fullfilled (for the longitudinal case, Eq. lfT^ . see Fig. Next, the measured Wi(k) 
functions were fitted with ^-independent effective masses. Finally their values were compared to the 
tree-level estimates: 

Ml^,f = Ml,,, = e'7, Ml,f = p. (20) 

Figj21 displays the measured equations of state of the three types of fields. It is obvious that the 
longitudinal and the transverse gauge fields show degenerate behaviour despite of the rather different 
formal procedure of the determination of and wt- The fitted squared masses appear in Tabled 







Lagrangian mass squared 


equipartition mass squared 


Tkinetic{\^\) 


(^totai{\m\^) 


scalar 


trans & long 


scalar 


trans 


long 


0.0409 
0.0102 
0.0025 
0.0006 


1.91 
0.47 
0.12 
0.03 


1.80 ±0.05 
1.95 ±0.03 
1.99 ±0.01 
2.00 ±0.00 


0.93 ±0.02 
0.98 ±0.01 
1.00 ±0.00 
1.00 ±0.00 


2.16 ±0.02 
2.23 ±0.02 
2.26 ±0.03 
2.26 ±0.03 


0.95 ±0.01 
0.99 ±0.02 
1.00 ±0.02 
1.00 ±0.03 


0.92 ±0.01 
0.98 ±0.01 
0.99 ±0.01 
1.00 ±0.01 



Table 1: The comparison of the averaged mass terms of the Lagrangian (e.g. M^ip^, ipi = p, At, Al) 
with the thermodynamically fitted effective squared masses Mf^^j. 

It is remarkable that for the Higgs field the measured values of the squared mass are 10-15% 
higher than expected on the basis of ()20|1 . This is the effect of the ultraviolet fluctuations which 
modify at one loop its mass [12] the following way: 

,.2 2 / 2 A\ 0.226 , , 

Ml^iatuce =^'+ + ( ± 2 1 —r- (21) 

Indeed, the variation of the deviation from the tree level squared mass at very low energy densities 
was found to be inversely proportional to a?. The non-perturbatively fitted coefficient is much 
smaller, than the 1-loop estimate appearing above. The convergence of the measured Higgs mass 
to a well-defined value for T {) makes it possible to define the renormalised Higgs mass. If one 
chooses it to coincide with the curvature of the classical potential at its minimum (e.g. then 
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Figure 3: Correlation between the scalar field and the longitudinal and transversal gauge fields. On 
the left the binned averages of ()22|) are displayed as a function of time (measured in units |m|~^). 
The error bars show the actual amplitude of its fiuctuations. The average values of the correlation 
coefficients are shown in the right hand figure as a function of the energy density (measured in units 
|m|^). The variation is compatible with the vanishing of the correlation for vanishing energy density 
(its slope is: 1.04 ±0.04). 

with the same subtraction one can define also the temperature dependent mass. In this way a very 
good agreement of the renormalised squared masses with the predicted temperature dependence was 
found. 

The effective squared masses of the vector modes were found very close to My = e^p^ without 
observable lattice spacing dependence. 

3 The statistics of the elementary fields 

In this section we shall analyse the data in further details focusing on the statistical correlations of 
the field variables p, A^,, Ay. Our aim is to find the field-coordinates of the quasi-particles, emerging 
from the thermodynamical analysis of the previous section. There are two aspects of this search: 

• Can one consider the above variables mutually independent in statistical sense? 

• Do they follow Gaussian statistics expected for small amplitude, nearly free quasi-particle 
oscillations? 

As a signature of independence of two field variables we shall consider the factorisability of the spatial 
average of their quadratic products by testing if, for instance. 



. ^ p^(x,t)A^(x,t)-p^(x,t) A^(x,t) ^ ^ 
^ p2(x,t)A2.(x,t) 

is fuUfiUed (analoguous test is performed for A[Al, p]). For variables obeying Gaussian statistics the 
equality 

(A2,(x,t))'-3(A2,(x,t))' 
r[Ar] = ^ ^ = 0, (23) 

(A|^(x,t))^ 

is obeyed. 
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Figure 4: Left: Three examples of the non-equihbrium evolution of A[Aj^,p]. In two time-histories 
stationary nonzero values of ~ —0.6 occured for a restricted time interval. It signals the presence of 
a vortex-antivortex pair and can be easily distinguished from the fast and smooth relaxation to zero, 
when no vortices are generated. Right: The vortex-antivortex pair corresponding to non-zero A, is 
found by looking for the points where p < 0.3 

Clear quantitative evidence appears on the left hand side of Figl^l for the stationary (time- 
independent) independence of and p. On the other hand and p show non-zero correlation, 
which appears also constant in time. The test of Gaussianity shows that from very early times both 
r[AT] and V[p] vanish, while r[AL] violates the criterium concerning the Gaussian nature of its 
statistics. 

On the right hand side of Fig. Elwe give the dependence of A[A2,, p] on the average energy density 
of the system. (The corresponding temperatures, understood in the sense described in section 2, are 
listed in Table 1). The effect increases linearly with the energy density. Analogous behaviour is 
observed for Fp^] and A[n2,,p], where TIl is the canonical momentum field conjugate to A^. 
The correlation coefficient A(Ai,p) was computed also on lattices of the same physical size but 
smaller lattice spacing providing essentially the same result. The lack of the p — factorisation 
was demonstrated on a number of functions /(p) replacing p^, among them also (p — p)^. 

In systems with large values of the couplings this strong correlation could be natural to expect. 
In this respect it is more surprising that the transversal gauge field behaves as a perfect independent 
variable. In non-equilibrium time evolution, when in addition to the small amplitude white noise 
initial conditions applied to the k 7^ modes, the average Higgs field (k = 0) starts from the unstable 
maximum at p = we have observed frequently the generation of rather large nonzero values for 
A [At, p]. Some spectacular examples are presented in FiglH One observes a rather abrupt vanishing 
of these values after certain time elapses. It turns out that in case of the transverse fields nonzero 
quasi-stationary values of A[At,p] are very sensitive indicators of the presence of Nielsen- Olesen 
vortex-antivortex pairs [23 • A representative vortex-antivortex pair is displayed in the right hand 
figure of FigEJ The drop of the absolute value of A[At, p] occurs when the pair annihilates. Similar, 
but less pronounced drop is observed in the correlation of the Higgs field and the magnetic plaquette 
variable, which is strictly gauge invariant even for finite lattice spacing. It is clear that in the vortex 
solution presented in the transversal vector potential and the real Higgs field vary in a correlated 
way. The detailed statistics of the vortex production during spinodal instability will be discussed in 
a separate publication. Here we conclude as a byproduct from this analysis, that the vanishing of 
A[At,p] in equilibrium not only hints to the independence of the A^-statistics, but also represents 
evidence that no vortex lines are present in our equilibrium configurations. 
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vve nave mvesxigaxea aiso me aepenaence oi p\ ana n.[j\T,P\ on xne niggs-gauge mass 

ratio, that is on A/e^. The value of A was varied in the range (0.75, 75) keeping the energy density 
fixed and = 1. A[Al,p] hnearly increases from ~ 0.001 up to ~ 0.1. Specifically, no effect on 
this linear variation can be observed when the decay channel of the Higgs into two "photons" opens 
around A/e^ ~ 12. The time average of A[At, p] stays zero, only the size of its fiuctuations increases. 
We have checked that the correlative effects do not change for = 0.01. 

In concluding this section we elaborate on the consequences of the statistical independence of At 
and p on the thermodynamics. The Gaussian nature of the Fourier modes of At is expressed by the 
following formula: 

AT(k)A^(k') = lArm^^y . (24) 

Let us apply it to the statistical expectation of a weighted space-average of At(x). One can make 
the observation that both and pt are quadratic expressions of At, multiplied by some function of 
p. First one substitutes for At its Fourier expansion. When the independence of the At statistics 
from p is used, the average can be written in a factorized way. Finally, is used and provides a 
simple general formula: 



^ J d3x/[p(x,t)]A2(x,t) = rf3a;/[p(x,t)] ^e^(k-k')x X AT(k,t)A^(k',t) 

k,k' 

= (/[p(x,t)])y Y: |AT(k,t)P, (25) 

k 

where the bracket with lower index V was introduced to denote the spatial averaging. This represen- 
tation allows for the following identification of the corresponding spectral density, with no ambiguity 
left in its definition: 

(/[p]A2)v(k) = {flpi^Mv |A(k,t)p. (26) 

Direct consequences of this are the following expressions for ZtiPt-, where factorized averaging over 
At and p appears: 



e^(k) = 1 [|nT(k)|2 + (k^ + (e2p2)v)|AT(k) 
1 



p^(k) = - |nT(k)|2 + (k^ - (eV)y)|AT(k)|2 . (27) 

They are equivalently given through the expression of the dispersion relation and the equation of 
state of the transversal quasi-particles: 



4(k) = affi = k^ + ^T(k) = ^ = 1 (28) 

^ |AT(k)|2 ^ ^ ' eT(k) 3k2 + (eV)v 

(The dispersion relation emerges from the average kinetic-potential energy equality for each mode). 
These expressions reproduce perfectly the original transversal spectral densities and the EoS ex- 
tracted in the previous section. 

Similarly good representation is obtained for the thermodynamics of the Higgs-sector when the 
Gaussianity of p is imposed on the averaging. In case of A^ a very bad quality representation of the 
spectral equation of state is possible if independence and Gaussianity are forced upon the longitudinal 
part of the vector potential. The quality of the description can be characterized by the violation of 
the kinetic-potential energy equality (|T^ . which is displayed in FigjHl 

Since the characteristics of the longitudinal equation of state clearly possess a quasi-particle 
character, as was seen in section 2, our conclusion is that A^ cannot be the corresponding quasi- 
particle coordinate. Instead in the next section we search for an appropriate non-linear combination 
of At and p. 
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Figure 5: Correlators of two versions of the nonlinearly mapped "longitudinal" field as a function of 
the parameter a. On the left the results of the mapping on the right the mapping (j21) are shown. 



4 Search for the longitudinal quasi-particle 

We search for nonlinear combinations of the longitudinal vector potential and p for which in (j^^ 
and (^Hj) equalities can be verified. This should support the general view of almost free quasi-particle 
composition of finite temperature media. It should be checked whether for this new combination the 
equipartition and the spectral equation of state can be interpreted consistently under the assumption 
of the factorisation of their statistics. It will be also checked to what extent it follows Gaussian 
statistics. Eventually we show that the quality of the equation of state based on this independent 
new variable is almost as good as the result of Section 2. We shall work out the example of (0), and 
only comment on the performance of the alternative proposition (j21). 

The coefficient a in (0) could be estimated analytically by assuming that it is small. For a 
quantitative treatment there is no need to apply approximations. By tuning a continuously we have 
measured the normalised reduced correlation functions A(^L,p) and A(Pl,p). In Fig. these two 
curves are displayed as functions of a together with the quantities testing the Gaussian nature of the 
new variables. It is remarkable that both the new coordinate and its conjugate momentum become 
un correlated with the Higgs-field for the same a. Although they are not perfect Gaussian variables, 
the deviation from Gaussianity is minimal just for this a value, it is about 5 times smaller than at 
a = 0. The situation is very similar for the nonlinear mapping (jSJ as can be seen from the twin 
figure of Fig. |S1 It is important to notice that in the latter case a 7^ 1, which would be the naive 
expectation on the basis of the way the kinetic-potential equipartition is fulfilled. Further in the 
analysis we use the optimal a values found as the average over 40 independent configurations. It 
has to be emphasised that no unique ao^t exists for single configurations, the optimal values found 
from the four different criteria fiuctuate somewhat. It is the ensemble average which leads to the 
remarkable coincidence displayed in Fig. |31 

Now one has to reexpress the partial energy densities and pressures in terms of the new variables. 
One should not forget that as a consequence of (H)) also the canonical momentum conjugate to p 
changes: 

^^^Tfb^'-'^'^'^' '''' 

Assuming factorisation of the statistical averaging over p. Hp and Al-,Vl one obtains the following 
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Figure 6: The quality of the mode- by-mode equahty of the longitudinal kinetic and potential energies 
when i) the original A^^, ii) the nonlinearly constructed variable A^^ are proposed as statistically 
independent fields. For comparison the perfect fullfillment of (fTK|) is also presented . The region 
k > 2\m\ is shown, where the equipartition is very well fullfilled for the other degrees of freedom. 

expressions for the averaged energy densities: 



+ , (30) 



2(1 + ap2 

It is worth now to check if the equality between the first two and the last term of Eq. ()30|) is fullfilled 
in the present factorised form. Fig. |B1 demonstrates no more than 1% deviation from the almost 
perfect equality found in Section 2. 

The corresponding expressions for the partial pressures are as follows: 



l («Vp^)^ p2 1 (1 + V )^: 



It is worth to write explicitly the spectral equation of state resulting for the new effective field 
variables, when the improved but still approximate equipartition is exploited: 



Wnew{k) = - ====. (32) 

3(l + V)2+(£5£!)! + k2(iW)! 
11 
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Figure 7: Sperctral equation of state with (continuous line) as compared with the curve 
("theory", dotted hne) where the p dependence of the coefficients is averaged numerically. The 
equation of state obtained in Section 2 ("without factorisation", dashed line) is also shown for 
comparison. 

The correction in the energy density and the pressure of p turns out to be negligible. In Fig. [3 
we plot the "longitudinal" equation of state using (jHOl) and (p?T|l . by measuring the corresponding 
averages directly, and the "theoretical" curve ()32|1 . A rather remarkable agreement of the curves is 
seen which works the best for low and moderate values of k. The systematic deviation experienced 
in the high k region would require smaller mass, which however would lead to large discrepancy in 
the low k region. Similar features are observed when the equation of state is constructed with help 
of the mapping Q. The curve corresponding to (jH^ agrees very well with the longitudinal equation 
of state appearing in Fig. El 

5 Conclusions 

In the present investigation we searched for the quasi-particles of the classical Abelian Higgs model 
at finite (low) energy density. The validity of a quasi-particle representation of the thermodynamics 
is essential for cosmological applications. The intuitively expected statistical independence and 
Gaussian behaviour of the transversal gauge field and of the Higgs excitations in the unitary gauge 
was confirmed. 

The longitudinal component of the vector field is not independent statistically at finite energy 
density from the Higgs oscillations. This correlation and the deviation from the Gaussian statistics 
grows with the temperature and also with the A/e^ ratio. Still, its thermodynamical characteristics 
represented in form of a spectral equation of state were shown to be degenerate with the features of 
the transverse vector modes. 

The first step of the search for an independent Gaussian pair of canonically conjugate field 
variables describing the non-transversal gauge fluctuations was put forward in this paper. We are 
confident that a fourth quasi-particle degree of freedom completes the thermodynamical character- 
isation, though there is no unique prescription for its systematical construction. One possibility is 
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to aoanaon lae laea or a locai mapping m x-space ana try xo use a K-aepenaenx mixing parameter: 
a(k). 

The present experience is very useful also for the investigation of the real time evolution of the 
field excitations in an infiaton+Gauge+Higgs system starting from an unstable (symmetric) initial 
state and ending in the broken symmetry phase. We are interested in the rate of excitation of the 
different field degrees of freedom immediately after the exit from infiation. Ah^cady at present it is 
clear that vortex-antivortex pair production will be rather frequent and can be studied also with help 
of the correlation coefficients introduced in the present study It will be interesting to see whether the 
energy transfer from the inflaton to the quasi-particles related to the transverse and the longitudinal 
vector fields is symmetrical. These questions are under current investigation. 

Appendix 

The thermal equilibrium of the system was realised in our real time simulations by long temporal 
evolution. It started from a noisy initial state having finite energy density. In order to implement 
the equations of motion in the temporal axial gauge, a spatial lattice was introduced. On this lattice 
the gauge field and the covariant (forward, F) derivatives of the scalar field are approximated with 
help of link variables: 

[/.(x) = e-^"l"*l^*(^\ D[^{x) = {Ui{x)^{x + tti) - ^{x)) . (33) 

The compact nature of C/'s has no effect at very low energy densities studied in the persent inves- 
tigation. The solution of the equations of motion was started with all conjugate momenta set to 
zero. Under this condition the Gauss constraint is trivially fulfilled in the temporal axial gauge and 
it is preserved during the evolution. The Fourier modes of the canonical coordinate fields were given 
an amplitude with random phase, providing for each mode equal potential energy. The simulations 
were done mostly on N — 32 three-dimensional lattices with a lattice constant a — 0.35|m|~^. In 
order to test the possible lattice constant dependence of the effects the system was solved also on 
= 64, a = 0.174|m|^^ lattice. Notable lattice spacing dependence was observed only for the Higgs 
mass (see the main discussion). 

The vector potential and its conjugate momenta are computed on an isotropic lattice (a^ = a) as 

A(x) = -^ImC/,(x), n, = -^Im (Ui{x + a5t)U-\x)) , (34) 

where the discrete time-step we use is = aSt with 6t ~ 0.1. These fields were decomposed into 
tranverse and longitudinal components, we need in the formulae for the decomposition of the energy 
and pressure. The magnetic field strength is measured from the plaquette variables 

Uij = Ur\x)U-\x + aj)Uj{x + ai)Ui{x), 
ImUij = -ea^eijkBk. (35) 

In the moments of measurements the actual configuration was transformed to unitary gauge. 
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